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The low energy quasiparticle dispersion of various narrow gap and gapless semiconductors are 
respectively described by three dimensional massive and massless Dirac fermions. The three dimen- 
sional Dirac spinor structure admits a time-reversal invariant, odd parity and Lorentz pseudoscalar 
topological superconducting state. Here we derive the effective field theory of this topological paired 
state for massless Dirac fermions in the presence of a fluctuating Zeeman term, which appears as 
a chiral gauge field. The effective theory consists of a mixed electromagnetic and chiral anomaly 
term in the bulk, and a combination of pure and mixed anomalies for the surface. In this paper 
we demonstrate the existence of fermion zero modes in the dilute vortex phase under generic condi- 
tions. Guided by the existence of the zero modes and its intimate connection with the anomaly, we 
propose an effective topological field theory in the presence of Dirac mass. We briefly discuss the 
experimental consequences of the effective field theory and the zero modes for the low temperature 
unconventional superconducting states of Cu a; Bi2Se3 and Sni-^In^Te. 



Introduction: Recently there has been a surge of the- 
oretical and experimental interest in the time-reversal 
invariant topological states of matter in three spatial 
dimensions |lH3[- In contrast to the abundance of mate- 
rials in which a three dimensional Z 2 strong topological 
insulator phase has been realized, there are only a few 
candidates for three dimensional, time-reversal invariant 
topological supcrfmid (TSF) and superconducting (TSC) 
states. Apart from the 3 He — B phase, which is a charge 
neutral TSF 0, 0] , there is only one charged, oddpar- 
ity TSC state, which has been proposed in Ref. |5| for 
Cu K Bi 2 Se3 [6|-l8|. The odd parity TSC state may also oc- 
cur in other narrow gap semi-conductors and semi-metals 
with a Dirac like quasiparticle dispersion in the normal 
state e.g, Su^^In^Te [9(. Both of these paired states 
belong to class DHI in the periodic table for topological 
insulators and superconductors [lot Hit . The BdG Hamil- 
tonians for uniform topological paired states in class DHI 
possess a nontrivial Z invar iant (TV) 0,1, Hlj . At an in- 
terface between two bulk states with different TV, the 
jump A7V governs the number of two component, two 
dimensional, massless Majorana fermions at the surface. 
The vortex physics of topological paired states is also 
supposed to be rich, due to the presence of gapless one 
dimensional modes along the vortex cores 

The recent point contact spectroscopy measurements 
on Cu2:Bi2Se3 I2-14| and Sni-^In^TeQ, have revealed a 
zero bias conductance peak (ZBCP), which is consistent 
with the existence of gapless Majorana surface states. 
However surface Majorana fermions may also emerge for 
other possible odd parity superconducting states, which 
are gapless in the bulk (such states are not described 
by the Z invariant of class Dili) JT2]. The magnetiza- 
tion measurements on Cu a Bi2Se3 also point towards an 
unconventional paired state jl5fll7| . On the other hand 
absence of ZBCP in a recent tunneling spectroscopy mea- 
surement has been interpreted in terms of a fully gapped, 
topologically trivial paired state [3] • Therefore more ex- 



periments arc needed to identify the pairing symmetry. 

So far, most of the theoretical studies have focused 
on the nature of the surface states in the presence of 
bulk topological pairing U 23| . The thermal Hall effect 



and its connection to gravitational chiral anomaly have 
been proposed as the appropriate topological response 
for the detection of the surface Majorana fermions (since 
charge is not a conserved quantity) in Refs. I24M26I On 
the other hand there are only a few theoretical studies 



27|,|28j of the bulk electrodynamic properties, which can 
also provide valuable insight into the nature of the under- 
lying paired state. In Ref. [27| a BF theory augmented by 
the fermion zero modes have been proposed. In Ref. |28| 
employing a higher dimensional (4+l)-d regularization, 
an axion electrodynamics action have been found as the 
effective field theory for (3+l)-d TSC state. However 
the importance of Zeeman coupling and the interplay be- 
tween charge and spin currents have not been considered 
in these papers. In this work we first consider the sim- 
pler case of massless Dirac fermions, and treat the orbital 
and Zeeman coupling on the same footing. By employing 
Fujikawa's chiral rotation technique p9j we derive the ef- 
fective field theory for the odd parity TSC state, which 
shows the presence of a mixed electromagnetic and chi- 
ral anomaly in the bulk, which becomes operative only 
in the vortex phase. We also find a combination of axial 
anomaly (pure electromagnetic + pure chiral anomalies) 
and mixed anomaly for the surface states. 

Normal state Hamiltonian: The low energy, long wave- 
length quasi-particle spectrum of several narrow gap and 
gapless semiconductors are respectively described by four 
component massive and massless Dirac equations in three 
spatial dimensions. In the presence of electromagnetic 
field, the text-book Dirac Hamiltonian is 

H N = J (fx^i-i^dj + e^Vf, + m)yb, (1) 

where Fermi velocity v has been set to unity, m is the 
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Dirac mass (band gap), e is the electron's charge, and 
V M is the electromagnetic vector potential. The anti- 
commuting 7 matrices satisfy {7^, 7"} = 2p MI/ , where 
guv _ ^ — 1 7 — 1 7 —1) is the metric tensor and %p = ^7°. 
Throughout the paper Greek and Latin indices will re- 
spectively correspond to space-time and spatial com- 
ponents of a vector. For simplicity we are ignoring 
the spatial anisotropy (direction dependent Fermi ve- 
locities) inherited from the crystal symmetry. If we 
choose the four component spinor ip T = (c^ ,cf , ojT , cj), 
where cf respectively correspond to the annihilation 
operators for parity even and odd states, with spin 
projections s 4,, the gamma matrices follow the 
so-called Dirac representation: 7 = n 3 <E> a , 7 J = 
in 2 ® cr- 5 , and 7 5 = i7°7 1 7 2 7 3 = n 1 ® a , where 
{7 5 ,7 M } = 0. The Hamiltonian H^ is invariant with 
respect to the parity (V), the time reversal (T) and 
the charge conjugation (C) transformations, respectively 
defined by Vij^t^V- 1 = 70^^,-x), Tipit,*)^ 1 = 
— 7i73V'(— t, x), and Cip(x)C" 1 = —i^2^*{x)- In the ab- 
sence of gauge field, the massless Dirac Hamiltonian is 
invariant under the global chiral transformation ip — } 
e «(0/ 2 )75^ ^ _^ ^ e %{e/2)y 5 _ rp^g S y mme t r y i s broken by 

the Dirac mass at the classical level, and by the gauge 
field through quantum mechanical effects known as chiral 
anomaly. 

Unlike the true relativistic field theory of electrons 
and positrons, the above Dirac description of the quasi- 
particles is just an emergent low energy theory. For elec- 
trons and positrons, the electromagnetic gauge field only 
appears in the Hamiltonian through covariant derivative, 
and this is also tied to the fact that gyromagnetic ratio 
for such fundamental particles is 2. For quasiparticles of 
a semi-conductor this is no longer true, and for this rea- 
son we need to consider the Zeeman coupling separately. 
When the high energy bands are integrated out in the 
presence of the electromagnetic field, the g factors (actu- 
ally a tensor) of two bands can be different [io| • For this 
reason the Zeeman coupling (again ignoring the spatial 
anisotropy) can be described by 



( 2 ) 



where 5+ ± respectively describe the g factors of the 
parity even and odd bands, and B = V x V is the mag- 
netic field strength. In the above equation we have ab- 
sorbed the Bohr magneton /is in the definition of g± . No- 
tice that the Zeeman term proportional to g+ appears as 
a vector potential for the chiral gauge field. We have in- 
troduced the matrices cr ji = {7^, 7 Z }/2i, and the g- term 
breaks chiral symmetry and charge conjugation symme- 
try. When the electromagnetic and the chiral gauge fields 
are simultaneously present, the chiral anomaly leads to 
the non-conservation of both chiral and electromagnetic 
currents (respectively defined by j^' 5 = ipj^^ip and 
j 1 * = V>7 m- 0)j which describes a novel interplay between 



orbital and spin currents. After integrating out the high 
energy bands we can also obtain different polarizabilitics 
for two bands and this effect can be captured through 
a term proportional to a 0: >Ej, where er - 7 = {7 ,7 J } /2i, 
and this term also breaks chiral symmetry. If we consider 
the simple problem of massless fermions and turn off all 
other chiral symmetry breaking bilinears (e.g., g- = 0), 
we obtain 
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where .F^/s an( i -^aS respectively correspond to the flux 



strength tensors of electromagnetic and axial gauge field. 
Notice that the non-conservation of electromagnetic cur- 
rent is captured by a mixed anomaly term in Eq. 01 and 
rather the modified electromagnetic current j a — j a — 
e 47T i •A-p^Yx wm be conserved. For a general chiral gauge 
field the modification corresponds to non-dissipative chi- 
ral magnetic current and anomalous Hall current (3li |32| . 
Since our chiral vector potential is proportional to the 
physical magnetic field, Ao = (temporal gauge condi- 
tion) . A "chiral electric field" will appear only for a time 
dependent physical magnetic field (B), and a spatially 
varying B can only induce a "chiral magnetic field" . We 
also note that a constant magnetic field will correspond 
to a constant chiral gauge potential, which corresponds 
to momentum space separation of left and right handed 
fermions in a Weyl semi- metal 31, 32|. In the following 
sections we will show that a mixed anomaly can occur 
even in the superconducting phase, but its effect can only 
be observed in the mixed phase. 

Pairing symmetries : In the presence of superconduc- 
tivity, we choose a Nambu spinor ty T = [tp^ , ip T j 5 C), and 
the 8x8 pairing Hamiltonian operator can be compactly 
written as 



H 



h N>1 - fit 

At 



hN,2 + A f l 



(5) 



where /ijv,i = 7 (— ij J, Dj + m + eAo7°), Hn,2 = 
7 (i7 J T>* — m — eAo7°) the covariant derivative T> j = 
(dj + ieAj) and Aj is the electromagnetic vector poten- 
tial. The pairing matrix A can be decomposed in terms 
of sixteen 4x4 matrices, which complete the Clifford al- 
gebra [33|]. The explicit form of the generalized pairing 
matrix is described by 

A = {A s t + + AV + AJ 7 " 7 S + A%a^}, (6) 

Under the Lorentz transformations (FT), the pairing bi- 
linears proportional to A 8 , A p , AJ^, A", A* respectively 
transform as scalar, pseudoscalar, vector, axial vector 
and antisymmetric tensor. For local pairing (intra-unit 
cell and momentum independent), the Pauli exclusion 
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principle only allows A s , A p and A^. The transformation erations, chiral transformation, and corresponding quasi- 



of the allowed pairings under LT, discrete symmetry op- 



particle spectra (for uniform paired states) are shown in 
TABLE I. 



TABLE I. Transformation properties of the allowed local pairing channels under the discrete symmetry operations V, T, chiral 
transformation U c , and the corresponding quasi-particle spectrum within the mean-field approximation. The distinct branches 
of quasi-particle spectrum are denoted by a = ±1, and each branch is two-fold degenerate. In Ref. [H the symmetries of 
possible local pairings for Cu^P^Ses, have been classified according to the representations of point group D^d- The relation 
between the decomposition in Our (A s l, Aq7°), A p i-y 5 , A37 3 , (A17 1 , A27 2 ) respectively correspond to A\ g , A% u , A 2u and E u 
representations of D-jd- 
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We also display the corresponding classification accord- 
ing to various representations of D^d point group, which 
was established in Ref. 0. 

Notice that only topologically trivial pairing A s is 
fully gapped for arbitrary values of the band parame- 
ters. In contrast the topological odd parity pairing A p , 



goes through a phase transition at \i + A 



(via 



a Dirac point) to a topologically trivial state. Only 
when /i 2 + A 2 > m 2 we have a nontrivial Z invariant, 
N = sgn(A p ) [33|. Consequently the Majorana surface 
states can only be found, when this condition is satisfied. 
For Cu x Bi2Sc3 and Sni-^In^Te, the chemical potential 
lies in the conduction band, this condition is always sat- 
isfied. The vector pairing A^ are not fully gapped. The 
space like vectors fi = 1,2,3 can have point nodes when 
/j 2 + A 2 > to 2 . For example if we consider j — 3, the 
nodes will be at k z = ±i//i 2 + A 2 — to 2 . The time like 
vector pairing exhibits gapless excitations on a three di- 
mensional surface in the momentum space. In the fol- 
lowing sections we will only focus on the fully gapped A s 
and AP. 

Effective action via chiral rotation: The derivation of 
the effective action in the presence of different chiral sym- 
metry breaking perturbations in the particle-hole chan- 
nel is a challenging task. However we can gain valu- 
able insight by first considering the pairing of mass- 
lcss Dirac fcrmions in the absence of any chiral sym- 
metry breaking perturbation in the particle-hole chan- 



nel. In the chiral basis of the four component spinor 



= (c: 



R „R 



), the massless Dirac Hamiltonian is 



block diagonal. The annihilation operators correspond- 
ing to the right (R) and the left (L) chiralities are re- 
spectively the symmetric and the antisymmetric combi- 
nations of the annihilation operators of the even and the 
odd parity bands cf^ L = (c+ ± c~)/ y/2. In this basis the 
trivial and topological pairings are respectively described 
by A s (RUa 2 R* + LHa 2 L*) and AP(RUa 2 R* + LUa 2 L*), 
where we have defined the two-component spinors R T = 
{cf,cf) and L T = (c£,cf). Therefore A s / p do not not 
mix R and L sectors and we can consider the pairing 
problem separately in R and L subspaces. The action 
in the fermion sector can be written as S = Sr + Si, 
the action in each subspace takes the following compact 
form, 



S a = / atx * a 



-A. 



*a (7) 



where a — R/L denote the chiral subspaces. For 
the chiral subspaces we have defined two four compo- 
nent Nambu spinors = {R,ia 2 R*) T and f L = 
(L,ia 2 L*) T , and two sets of 4 x 4 gamma matrices 7^. 
The explicit forms of these gamma matrices are given 



by 



1r = 1l= V 3 



a", % = -% = -ii] 2 <E> cr J , 
(7°. The right and left sectors now have 
separate gauge fields A^^ R / L = eV^ ± g+A v . Also notice 
that we have allowed two distinct phase angles Or/l- 
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In this decoupled chiral Nambu basis, the pairing term 
appears as a mixture of scalar and pscudoscalar Dirac 
masses, albeit in the presence of a new chiral gauge field 
in the Nambu space. This stems from the fact that (J2R* 
and 01L* respectively transform as left and right chiral 
spinors. Therefore in each chiral Nambu subspace, the 
phase of the superconductor appears as the axion angle 
and the gauge field appears as an axial gauge field. For 
trivial pairing Or = Ol = is the phase corresponding 
to electromagnetic U(l) symmetry. For the topological 
pairing Or — 9, 6l — ± 7T, where 6 again corresponds 
to the electromagnetic U(l) phase. 

For deriving an effective action of superconducting or- 
der parameter by integrating out the fermions, it is a 
common practice to first absorb the half of the supercon- 
ducting phase onto the fermion operators, which modi- 
fies the gauge field — > — 8^9/2. In our problem 
such a phase transformation ^f a — > exp {—i/20 a 7^) ^ a 
constitutes a chiral transformation and a modification 
of the chiral gauge field. In this process we should en- 
counter field theory anomalies, which leads to the topo- 
logical terms in the effective action. Following Fujikawa 
we perform the chiral rotation at infinitesimal steps 
v & a — > exp (— i/2 ds 8 a 7^) f n , and integrate over s from 
to I. The topological contribution arises through the 
Jacobians of the chiral transformation, and are given by 



J a = exp 



( I D 2 

~5 / ^a,s 



ds lim Tr < 7^ exp 



M—too 



A/ 2 



where the expression for the Dirac kernel is 



(8) 



D a =i^[d^ + \ A^ a 



^ _ Ae «{l-<is)««7S 



(9) 



The topological terms in this procedure is proportional 
to Tr[^a^o-? x }T^, a T p x, a . Now notice that 



tR = iiRiRiRiR, fL^^ihiiiii (10) 



and this difference leads to opposite signs for the topo- 
logical terms obtained from the right and the left sec- 
tors. After performing a standard momentum integral 
and combining the contributions from both sectors we 
obtain 
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(11) 



which is the main result of this paper. 

We first consider topologically trivial superconductor, 
and using 6r = Ol = in Eq. [TT] we obtain, 



S, 



top,s 



e p,vp\ 

16tt 2 



eg + J d^xBTl v Tf x 



(12) 



On the other hand for topological odd parity supercon- 
ductor, we substitute Or — 0l = if and find 



S. 



top,p 



d 4 x 



e p,vp\ 

16tt : 
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L J pv J pX 



9\^px + ^g + j; v T l 



P x 



(13) 



Notice that for massless Dirac fermions, both supercon- 
ductors have a mixed anomaly present in the bulk. The 
axion angle 9 for the mixed anomaly is a dynamic phase 
variable, and will modify the Maxwell's equations in the 
bulk. This represents the interplay between charge and 
spin currents. In addition to the bulk term, the odd par- 
ity pairing has an addition anomaly term with constant 
axion angle 7r, which is inoperative inside the bulk. This 
term only becomes important at the surface, and cor- 
responds to the existence of surface Majorana fermions. 
The BdG nature of the quasiparticles has contributed by 
a multiplicative factor of 1/2 to the conventional formula 



of anomaly. Also notice that mixed anomaly term does 
contribute to the surface action. If we do not consider 
the Zeeman term (chiral gauge field), there will be no 
bulk anomaly, and only surface anomaly term for the 
electromagnetic gauge field will be present for odd parity 
pairing. In that case our Eq. [T3]will simplify to 

S t o P , P = --^T- J d'x^x (14) 

which precisely agrees with the topological term derived 
in Ref. I28I using (4+l)-dimensional regularization. In 
the topological insulator where electromagnetic charge is 
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a conserved quantity, the existence of such an axion term 
in the electrodynamics (with twice the coefficient) corre- 
sponds to half-inte ger quantum Hall effect of the sur- 
face Dirac fermions [351 ] . However in the superconductor 
electromagnetic gauge symmetry is broken, and electro- 
magnetic charge is not a conserved quantity. Therefore, 
there will be no surface quantum Hall effect in the present 
problem. 

The conventional terms of the Landau-Ginzburg the- 
ory, describing the Meissner effects etc. is obtained by in- 
tegrating out the transformed fermions, assuming slowly 
varying chiral supercurrents (2A^ a — d^9 a ) and A(x). In 
the London limit we can neglect the fluctuations of the 
amplitude A(x), and obtain 



S, 



GL 



{d li e-2eV li f+^lAl 



GO 2 

(15) 



In the above equation we have also ignored the mismatch 
between Fermi velocity and speed of light. Now combin- 
ing Sql with S top we obtain the total effective action 
for the superconducting order parameter in the London 
limit. Now taking a derivative of the total action in the 
bulk with respect to V M , we find that the conventional su- 
percurrent 2ep(d fl 9 — 2eV tl ) is not conserved, when we are 
in the vortex/mixed phase. It turns out that topologi- 
cal superconductor has fermion zero modes in the vortex 
core, which provides an additional path for transport, 
and the mixed anomaly is precisely tied to this effect. 
Motivated by this we now proceed to explicit compu- 
tation of fermion zero modes in the dilute vortex limit. 
Fermion zero modes in the absence of Zceman coupling 
have been recently addressed by various authors. How- 
ever we will account for the general Zeeman coupling 
and its possible variation due to Meissner screening away 
from the vortex core in the following section. 

Vortex zero modes of generic Dirac Hamiltonian with 
odd parity pairing : We consider a static line vortex along 
the z direction, and solve the zero mode problem for the 
planar part. In the absence of Dirac mass, chemical po- 
tential, and the Zeeman couplings the planar problem 
corresponds to two copies of Jackiw-Rossi Hamiltonians, 
which belongs to class BDI and we obtain 2n number of 
Majorana zero modes, where n is the vorticity (36j . In 
the presence of the additional perturbations mentioned 
above, the problem belongs to class D and only for odd 
vorticity we can find 2 Majorana zero modes [371 ] . which 



38, 3 



is tied to a Z 2 index theorem 

Let us consider a particular profile of the electromag- 
netic gauge potential V</> = -^r if r < A, defining the 
core of the vortex and = l/2r outside (r > A), 
if one wishes to commit to the spherically symmetric 
gauge. Then the magnetic field (B) is finite and constant, 
B = 2^3- only inside the vortex core, while it vanishes 
outside. The superconducting order parameter vanishes 



smoothly inside the core as r 0, and satutares at Ao 
as r — > oo, otherwise arbitrary. Upon including the Dirac 
mass (mfc), Zeeman couplings (hi, h 2 ) one set of coupled 
differential equations of zero energy modes in the plane 
perpendicular to the applied field read as 

e~** (d z + V4,) AT/ + ia(r) A+ + iA^"^ = 0, 

ie* (d g - V ) A+ + b(r) Al + A r e"^ (A+) f = 0(16) 

at finite chemical potential (fj,), where d z = d r — -d^, 
d z = d* and a(r) = mu + fJ, + hi(r), b(r) = m k + h 2 (r) - 
(i, and in what follows we set mk = rn (constant). In 
terms of g-factors, which can in principle be different in 
two bands, one can write hi(r) = p,£>(r), where i = 
1,2. Therefore, the Zeeman coupling is only finite inside 
the vortex core, and vanishes outside. The other set of 
coupled differential equations for the zero mode is 

e-** (d z + V ) A+ - id(r) A" + iA r e~^ (a+) ^ = 0, 
(i)e l<i> (d- z - V ) A- - c(r) A+ + A r e"^ (k^ = 0(17) 

where c(r) = mk — hi(r) + fi, d(r) = — h%(r) — p. 
The remaining four equations are related to these by 
Hermitian conjugation. 

Two BdG-Majorana modes in Nambu-Dirac basis, de- 
#!(+£), , where = 



fined as ^> = 



v& T (fc) and *^(— k) = ^p(k), otherwise 
* T (x) 

take the form 
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(18) 
\ 



(19) 



where TZ(r) = exp (-if - /J" A r ,dr'). For r/2A < 1 



g(r) = ci h 



ry ab 



, f(r)=ci 




ry ab 



(20) 



where Ij, with j — 0, 1 are the modified Bcssel functions 
of order j. On the other hand, deep within the vortex 
core 



q(r) = ci h 



P( r ) = ci \j- Io 



led 



(21) 
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Outside the core of the vortex (r > A) 

g(r) = c 2 I 1/2 (r^/Pm + c 3 J_ 1/2 (ry/prj 



f(> 



C2I-1/2 [ry/M + csh/2 [ry/fa) (22) 



where (3 = m + fi, and r\ = \m — In the regime, 
g(r) = q(r) and f(r) = p(r). Two out of three arbitrary 
constants, for each of the BdG-Majorana solutions can 
be fixed by equating the value of the functions, and their 
first derivatives at r = A. The remaining constant can 
be determined from the normalization condition. The 
modified Bcsscl functions I+i /?jar) oc e Qr /ra, grow ex- 
ponential at large distances [4C 



4ll | . Therefore, for nor- 



malizable zero energy modes can only be found if the 
condition 



Aq > mf 



(23) 



is satisfied. It is worth noticing that the nomalizability 
of the solutions is independent of the Zeeman couplings. 

We now account for the 83 and A3 along the vortex 
line perturbativcly. This perturbation reads as rl3 = 
(r (—ids) + T3A3) <g> 7°7 3 and keeps the BdG-Majorana 
sub-space invariant. In this basis H.3 is purely off diago- 
nal, and takes the form 



U3 = -Ui(T 2 (i) (d 3 - iA 3 ) - v 2 a 1 {i) 



d 3 1 A 3 

V2 



(24) 

Two parameters (v±, v 2 ) of the one dimensional helical 
BdG-Majorana modes depend on the exact profile of the 
vortex, and other microscopic parameters (Dirac mass, 
chemical potential etc.). However, in the absence of any 
Zeeman coupling, i.e. h\ = h 2 = 0, one finds v\ = v 2 . 
The above Hamiltonain in Eq. [52] can be diagonalized in 
the basis of symmetric and anti-symmetric combinations 
of I*?) and \\V%), yielding 



n. 



-V1W3 (d 3 - iA 3 ) 



(25) 



Notice that gaguge field is coupling in conventional way 
(unlike Nambu space), which reflects the quasiparticle in- 
side the core describes normal state. Since we have both 
left and right movers (in one dimension) coupled to the 
electromagnetic field, we do not have any electromagnetic 
anomaly in one dimension. This justifies that there is 
no builk anomaly in the absence of Zeeman/chiral gauge 
field. In the presence of the Zeeman coupling we find the 
both electromagnetic and chiral gauge field couplings for 
the one dimensional modes, which leads to the one di- 
mensional version of the chiral anomaly. The presence 
of the bulk anomaly term is crucial to feed the current 
to one dimensional world via Callan-Harvey mechanism 
42|. 

In the absence of the chiral symmetry breaking per- 
turbations (Dirac mass, g_) in the particle hole channel, 



fermion zero modes are also present for the topologically 
trivial pairing A s . However these zero modes get re- 
moved by such chiral symmetry breaking perturbations 
[? ]. Therefore we believe the mixed anomaly term will 
be absent for trivial superconductor under generic condi- 
tions. In contrast the zero modes of the topological pair- 
ing are robust under generic conditions, and this suggests 
that mixed anomaly term will be present for topological 
pairing, even in the presence of generic perturbations. 
However the coefficient instead of being eg+ will have 
more involved model dependence. In contrast a Dirac 
mass in the absence of time reversal symmetry break- 
ing perturbations, will modify the surface term by a step 
function 0(/z 2 + A 2 — m 2 ), which is the condition for 
finding masslcss surface Majorana fermions. The pres- 
ence of time reversal symmetry breaking perturbations 
will induce a gap for the surface Majorana fermions. 

Experimental consequences : We have found a general 
effective action for the topological superconducting state, 
which suggests an intriguing interplay between charge 
and spin currents. The bulk anomaly terms (propor- 
tional to the dynamic phase variable) will have nontriv- 
ial effects on the electrodynamics equations in the vortex 
phase. However these will be appreciable only for time- 
dependent electromagnetic fields. The derived surface 
anomaly terms will also play important role in the elec- 
trodynamic properties of Josephson junction measure- 
ments. The presence of robust fermion zero modes in the 
dilute vortex limit suggests a T-linear specific heat in the 
mixed state. The anomaly terms in the effective action 
suggests that such behavior should be present through- 
out the mixed phase. More detailed derivation of the 
effective action under generic conditions, and a thorough 
analysis of the electrodynamics of the mixed state will be 
presented in a forthcoming publication (43j . 
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